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Abstract 

The quantum resonances (QRs) of the kicked particle are studied in a most 
general framework by also considering arbitrary periodic kicking potentials. 
It is shown that QR can arise, in general, for any rational value of the Bloch 
quasimomentum. This is illustrated in the case of the main QRs for arbitrary 
potentials. In this case, which is shown to be precisely described by the linear 
kicked rotor, exact formulas are derived for the diffusion coefficients deter- 
mining the asymptotic evolution of the average kinetic energy of either an 
incoherent mixture of plane waves or a general wave packet. The momentum 
probability distribution is exactly calculated and studied for a two-harmonic 
potential. It clearly exhibits new resonant values of the quasimomentum and 
it is robust under small deviations from QR. 
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I. INTRODUCTION 



There has been recently a considerable experimental [1] and theoretical [2,3] interest 
in new remarkable phenomena associated with the quantum resonances of the periodically 
kicked particle (KP), either in the presence or in the absence of gravity. In this paper, 
we report about new general results concerning the quantum resonances of the KP in the 
absence of gravity. These results extend some of previous work [2,3] in a significant way, 
shedding light on new basic aspects of the problem and connecting an important case of the 
system with a well-known model. We start with a brief summary of previous results [2,3], 
using notation in Ref. [3]. The quantum KP is described by the Hamiltonian 

H = ^ + kV(x)J2S(t'-tr), (1) 

t 

where (x, p) are the position and momentum of the particle, k is a nonintegrability param- 
eter, V(x) is a periodic potential, t takes all the integer values, t' is the continuous time, 
and r is the kicking period. The units are chosen so that the particle mass is 1, the Planck's 
constant h — 1, and the period of V(x) is 2n. In practice, V(x) has been always chosen 
as the standard potential V(x) = cos(x). The one-period evolution operator for (1), from 
t' = t + to i! = t + t + 0, is given by 

U = exp [-ikV(x)] exp (-irf/2) . (2) 

The translational invariance of (2) in x implies the conservation of a quasimomentum (3 
(0 < (3 < 1): The application of U on a Bloch function ^p(x) = ex-p(i/3x)ifj/3(x), where 
^p(x + 2n) = ippix), results in a Bloch function = exp(i/3x)^(x) associated with the 

same value of (3. Here ip'^x) is the 27r-periodic function ip'pix) = (jpipp{x), where 

= exp [-ikV(x)} exp [-it (p + f3) 2 /2] . (3) 

The restriction of the operator (3) to 27r-periodic functions i^p{x) allows one to interpret x as 
an angle 6 and p as an angular-momentum operator TV = —id/d6 with integer eigenvalues n. 
One can then view (3) as the one-period evolution operator for a "/5-kicked rotor" {(3-KB). 

2 



Now, an arbitrary KP wave packet ty(x) can be always expressed as a superposition of Bloch 



for integer "time" t, connecting the quantum dynamics of the KP with that of /3-KRs. 

For typical irrational values of rj (2ir), /3-KRs are expected to feature dynamical localiza- 
tion in the angular momentum n [4,5], implying a similar localization of the KP wave packet 
^(p) in the momentum p = n + j3 (n and (3 are, respectively, the integer and fractional parts 
oip). If r/(2n) is rational, the usual (j3 = 0) KR exhibits quantum resonance (QR), i.e., a 
ballistic (quadratic in time) growth of its kinetic energy [5]. QR in general /3-KRs appears to 
have been studied only in the case of integer rj (2n) ("main" QRs) with V(x) = cos(x) [2,3]. 
It was found [3] that QR arises in this case only for special "resonant" values of (3, finite 
in number. Thus, QR is exhibited by the KP only if is delta localized on the discrete 
set of momenta p = n + j3, with (5 in the finite resonant set. However, QR leaves a clear 
fingerprint in the evolution of either a general KP wave packet (5) or an incoherent mixture 
of plane waves: In both cases, which involve all values of (3, the average kinetic energy grows 
diffusively (linearly) in time [2,3]. This diffusive behavior is robust under small deviations e 
of t/(2tv) from integers, in the sense that it is still observed on time scales t oc |/ce| -1 / 2 [3]. 

In this paper, the results above concerning the QRs of the system (1) are extended by 
also considering arbitrary periodic potentials V(x). In Sec. II, we show that the general 
condition for QR in /3-KRs is just the rationality of both t/(2tt) and (3. Thus, a resonant 
value of (3 may be a general rational number in [0, 1) and one then has a (countable) infinity 
of such values. This is illustrated in Sec. Ill in the case of the main QRs for arbitrary V(x). 



functions, ty(x) = d/3exp(i/3x)ipp(x), where 




(4) 



ty(p) being the momentum representation of *ff(x). One then gets the basic relation 




(5) 
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We show that this case is precisely described by the well-known linear kicked rotor (linear 
KR) [9,10]. We focus on this case also in Sees. IV and V. In Sec. IV, we derive exact 
formulas for the diffusion coefficients determining the asymptotic evolution of the average 
kinetic energy of either an incoherent mixture of plane waves or a general KP wave packet. 
For a mixture uniformly distributed in (3, the formula essentially coincides with one given 
by Berry in the context of the linear KR [10]. Multi-harmonic potentials cause the diffusion 
coefficient for a general wave packet to be explicitly dependent on quantum correlations 
reflecting the localization features of the initial wave packet in momentum space. In Sec. V, 
the momentum probability distribution is exactly calculated and studied for a two-harmonic 
potential. It clearly exhibits new resonant values of (3 due to the second harmonics and it is 
found numerically to be robust under small deviations of r/(2n) from integers. A summary 
and conclusions are presented in Sec. VI. 

II. GENERAL QR CONDITIONS 

The basic origin of QR is a band quasienergy (QE) spectrum due to some translational 
invariance in phase space. For KR systems, this is the invariance of the evolution operator 
under translations T q = exp(—iq9) by q in the angular momentum N [5,6]; here q must 
be an integer since 9 is an angle (0 < 9 < 2n). The translational invariance of /3-KRs is 
expressed by [Up, T q ] = 0, where Up is given by Eq. (3) (x — > 9, p — > N). Using the fact 
that N has integer eigenvalues, one easily finds that [Up, T q ] = is satisfied only if 

= j - § mod(l), (7) 

where / and r are integers. Eq. (6) is the rationality condition for t/(2tt) while Eq. (7) is a 
formula for the general resonant values of (3. For definiteness and without loss of generality, 
we assume that I and q are positive. Let us now write / = glo and q = gqo, where Iq and 
g are coprime positive integers and g is the greatest common factor of (I, q); the value of 



r/(2n) = lo/qo will be kept fixed in what follows. It is then clear that (3 in Eq. (7) can take 
any rational value /3 r in [0, 1) since g can be always chosen so that r = (/3 r + gq /2)gl is 
integer. For given (3 = f3 T , we shall choose g as the smallest positive integer satisfying the 
latter requirement, so as to yield the minimal values of / = gl and q = gq . In general, 
g > 1, so that (/, q) are not coprime. For the usual KR {(3 — 0), g — 1 if l q is even 
and g — 2 if Z <?o is odd (compare with Ref. [6]). We denote (3 X by /3 rj9 , where the integer 
r = (/3 r + gq /2)gl labels all the different values of (3 r for given minimal g. 

The QE states (p for (3 = (3 r ^ can be chosen as simultaneous eigenstates of Up and T q : 
Upip = exp(—iu)(p, T q ip = exp(—iqa)ip, where uj is the QE and a is a "quasiangle" , varying 
in the "Brillouin zone" (BZ) < a < 2-n/q. One may view the Bloch function exp(i(3x)(f(x) 
as a state on the "quantum torus" < x < 2n, < p < q, with toral boundary conditions 
[7] specified by (a, (3). Using standard methods [5,6], it is easy to show from the eigenvalue 
equations above that at fixed a one has precisely q QE eigenvalues oj^a, (3), b — 0, ... , q — 1. 
Since q = gq is minimal, the BZ is maximal for the given value of (3 — f3 r>g . Then, as a is 
varied continuously in the BZ, the q eigenvalues form q QE bands. These bands are expected 
to be, typically, not all flat (with zero width); QR can then arise and (3 = (3 r ^ g is indeed a 
resonant value. In the nontypical case that all the bands are flat, (3 = (3 r , g is nonresonant: 
QR is replaced by a bounded quantum motion, the "quantum antiresonance" [8]. 

III. CASE OF MAIN QRs: CONNECTION WITH THE LINEAR KR 

Previous studies [2,3] have focused on the important case of the main QRs (r = 2nl Q , 
q = 1), assuming that V(6) = cos(0). It was found [3] that only l values of (3 are resonant 
(exhibit QR). They are given by Eq. (7) with r = 0, 1,..., Iq — 1 and g = 1 (i.e., 
/ = Iq and q — 1). In this section, we study the main QRs for arbitrary V{6) and we 
show that all rational values of (3 in [0, 1) are resonant if V(9) contains all the harmonics. 
The case of r = 2tt1 is the only one in which the term (N + (3) 2 in Eq. (3) (x — > 9, 
p — > N) can be replaced by the operator N + 2(3N + (3 2 , linear in N; this is because 
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exp(— inl n 2 ) = exp(— iirl n) for the integer eigenvalues n of N . Then, after omitting the 
nonrelevant constant phase factor exp(— iirlo(3 2 ), one can express (3) as follows: 

Up = exp -ikV(§) exp (-irpNj , (8) 

where 773 = tiIq(2(3 + 1). We identify Up in Eq. (8) as the one-period evolution operator for 
the well-known linear KR [9,10] with Hamiltonian H = TpN + W{9) 'EZ-oo S ( t ' ~ *)l the 
corresponding Schrodinger equation is exactly solvable for arbitrary potential V(9), 

V(9) = J2V m exp(- im 9). (9) 

m 

Assuming the quantum state of the linear KR to be initially (at t = 0) an angular-momentum 
state, ippfl(9) = exp(ino9)/y/27i, the expectation value of the kinetic energy at time t is given 
by the exact expression [10,11]: 

If (5 is a typical irrational number, so that also Tp/(2ir) = l ((3 + 1/2) is such, the QE 
spectrum of the linear KR is pure point [9] and E n0t p(t) in Eq. (10) is a bounded, quasiperi- 
odic function of time [10]. Consider now a general rational value of f3 — (3 r , g , given by Eq. 
(7) with I = gl , q = gq = g, and g minimal (see Sec. II). It is easy to show that f3 r>g 
can be written, up to a fixed integer shift in r by gla/2 if both g and gla/2 are even, as 
Pr,g — r /(glo) ~ 1/2 mod(l), where r and g are coprime; the corresponding value of Tp/(2ir) 
is r/g, up to some additive integer. Then, from work [9], the QE spectrum consists of g 
bands which are all nonflat (i.e., f3 = (3 r>g is resonant) if there exists at least one integer 
j 7^ such that the Fourier coefficient Vj g in Eq. (9) is nonzero; otherwise, all the g bands 
are flat {(3 = (3 rj9 is nonresonant, with QR replaced by quantum antiresonance [8]). In fact, 
if Vjg 7^ for some j ^ 0, one finds from Eq. (10) (with Tp = 27ir/g) a ballistic behavior 
for large t, E n0y p{t) w St 2 /2, where S = k 2 g 2 Y Jj f\V jg \ 2 [10]. Thus, if V{9) contains all the 
harmonics (V m ^ for all m), all rational values f3 r>g of f3 are resonant. 
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As a simple example, let V(6) = cos(9) + 7cos(2#). The only nonzero coefficients V m are 
V±i = 1/2 and V± 2 = 7/2, so that QR arises only for g — 1, 2. The resonant (3 values are 
Pr, g = r /(glo) — 1/2 mod(l), where r = 0, 1, . . . , 1$ — 1 for g = 1 and r = 1, 3, . . . , 2/ — 1 
for g = 2. The values for g — 1 are just the known ones for V(9) = cos(#) [3] while those 
for g = 2 are new ones, due entirely to the second harmonics. See, however, note [12]. 

IV. ASYMPTOTIC DIFFUSION OF AVERAGE KINETIC ENERGY 
A. Incoherent Mixture of Plane Waves 

Let us assume, as in experimental situations [1], that the initial KP state is an incoherent 
mixture of plane waves exp(ipx) with momentum distribution f{p) sufficiently localized in 
p. By decomposing p into its integer and fractional parts, p = n + (3, the average kinetic 
energy of this mixture at time t can be expressed as E{t) = d(3 Yin f( n + P)^'n,p^)'i nere 
E' n f3 {t) is the expectation value of the kinetic energy in the state evolving from a plane wave. 
In the case of r = 27rZ , on which we shall focus, E' n Jt) is given by the right-hand side of 
Eq. (10) with n replaced by p = n + (3; this can be easily seen from Eqs. (18)-(21) in Ref. 
[10]. As in Ref. [3], we define fo((3) = J2 n f( n + P) an d use the asymptotic (large t) relation 

f 1 ^, m sin 2 [rcmlp ((3 + l/2)t] t 1 ^ 
Jo sm [nmlo {(3 + 1/2)] \m\l 

where f3 r ^ m = r/(\m\l ) — 1/2 mod(l). We then find that E(t) behaves diffusively for large 
t, E(t) ~ D t, where the diffusion coefficient D is given by 

,2 mlo— 1 

D = T Y,m\V m \ 2 fo (Pr,m) ■ (11) 

u m>0 r=0 

The simple mixture with f{n + (5) — 5 n>no is uniformly distributed in (3, fo((3) = 1. In this 
case, E(t) = J* d(3E' nQ p{t), leading to an exact equality for all times t: E(t) = E(0) + D t 
with D = A; 2 /(47r) jQ W [dV(9)/d9] 2 d9. The latter result was essentially obtained by Berry 
[10] as the kinetic energy of a linear KR whose value of r is completely unknown. 
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B. General Wave Packet 



Next, we consider a general KP wave packet (5), a "coherent mixture" of Bloch waves 
exhibiting all values of (3. The expectation value of the kinetic energy in ^(x) = (j t ^(x) is 



dd 



dd 



(12) 



where the last relation, with ipp,t{9) = UgippiO), holds for large t provided that (E) t is 
unbounded as t — > oo [2]. In fact, we now show that (.E) t exhibits an asymptotic diffusive 
behavior for r = 2nl . In this case, we easily obtain from Rels. (8) and (9) that 

M(9) = t)J^(0) = exp [-ifcy At (0)] MO - trp), (13) 

where 

%M = f.V(i- «,) = £ ^ 4 n n ( (^/ / 2) e " , '" l °" ( '"''" /21 - (14) 

As shown in Ref. [2] (Appendix A) for V (9) = cos(6 l ), with straightforward extension to 
arbitrary V (9), the expression (13) implies that for large t the dominant contribution of 
\di/j l3 ,t(9)/d9\ 2 to (E) t in Eq. (12) is k 2 |^(0 - *r^)trt^ )t (0)/d0| 2 . This contribution can be 
calculated using Eq. (14) and a relation following from Eq. (4): 

\M8)f = ^^Cp(m)exp(im0), (15) 

m 

where Cp{m) are correlations in momentum space, 

Cp{m) =^2^{m + n + I3)^!*{n + 13). (16) 

n 

From Eqs. (14) and (15), one can write the Fourier expansion of k 2 \^p(9 — tr^dVp^) /d9\ 2 . 
After inserting this expansion in Eq. (12), we get 

( E )t ~ % E mm ^C5(m, m'; t), (17) 
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where 



B(m to'- t)= C d Q Mmrpt/2) sin(mV/2) _ , K _ m)(t+1)W2 

1 ' ' tj 7 ^ sin(mr /3 /2) sin(m'r /3 /2) M J ' 1 J 

We show in the Appendix that the asymptotic behavior of the quantity (18) for t — > oo is 
given by 

g(m,m')l Q -l 

B(m, to'; t) ~ — (|m| + |m'| - |to - m'|) E Cp rg (m-m'), (19) 

^ r=0 

where (7 = <?(to, to/) is the greatest common factor of (|to|, |to'|) and /3 rj9 = r/(gl ) — 1/2 
mod(l), r = 0, . . . , g/ — 1. Rels. (17) and (19) imply a diffusive behavior of {E) t for large 
t, (E) t ~ Dt. The diffusion coefficient D can be expressed, after some algebra, as the sum 
of two terms: 

, 2 mio-l 

D = Dl + Du = T J2m\V m \ 2 C Pr,JO) 

m>0 r=0 
„, 2 00 00 g(m,m')/ -l 

+— E E mRe (KJO) E Re[C^ s (m-m')]. (20) 

m=l m'=m+l r=0 

For a potential containing harmonics V m of sufficiently high order to, the term Z?n [on the 
second line of Eq. (20)] can generally lead to a sensitive dependence of D on the correlations 
Cp(m), which reflect the profile of the initial wave packet in momentum space by Eq. (16). 
This term vanishes in some cases, e.g., for V(6) = cos(#) and/or a uniform probability 
distribution \ijjp{9)^ in Eq. (15). One is then left with only the first term (Di), which 

is completely analogous to the diffusion coefficient (11) for the incoherent mixture; this is 

~ 2 

because Cp(0) = ty(n + /3) from Eq. (16) and this is analogous to fo((3) = f( nJ <~P)- 
V. MOMENTUM PROBABILITY DISTRIBUTION 

The momentum probability distribution (MPD) for a KP wave packet is given by 

~ 2 „ 

P(p, t) = tyfip) , where ^tip) is the momentum representation of the wave packet at 
time t. At fixed /3, P(n + fl, t) is the angular-momentum (n) distribution for a /3-KR. Under 
QR conditions and for resonant /3, the motion of the /3-KR is ballistic in n, so that the width 
An(t) of P(n + (3, t) increases much faster than that for a nonresonant (3 [for most values of 



P, An(t) is expected to be essentially bounded due to dynamical localization]. Thus, if p and 
t are sufficiently large, P(p, t) is almost zero except of narrow peaks around p = n + (3 r;9 , for 
all resonant values (3 T , g of (3. The diffusion coefficient (20) is the average of P(p, t)p 2 /(2t) 
over p in the limit of t — > oo; the only nonzero contributions to this average come from the 
resonant peaks. The contribution of all the peaks with fixed (3 r>g = (5 r ^ g i is precisely the 
contribution of all the terms with (3 r , m = j3 r ^ g ' and/or (3 r , g = (3 r i, g i in formula (20). 

We now study the MPD for r = 2irl and V(x) = cos(x) + 7cos(2x) (see also end of 
Sec. 111). From Rel. (4), \l/ t (n + (3) are just the Fourier coefficients of \f2m^)p^{9) an d, 
for r = 27r/ , ipp,t(0) can be calculated using Rel. (13). We assume an initial wave packet 
ty(p) satisfying \£r(p) = 1 for < p < 1 and 4f(p) = otherwise. This corresponds, 
by Rel. (4), to a uniform ipp(9), ipp(9) = (2tt)~ 1/2 for all (3. Rel. (13) then implies that 
V2nipp,t(0) — ex P [—ikVp tt (0)] . Since only terms with \m\ — 1, 2 appear in the sum (14), the 
Fourier coefficients of exp [— ikVp :t (9)] can be easily expressed, essentially, as a convolution 
of ordinary Bessel functions J n (-) in the index n. We finally obtain the exact expression 



where kj(/3, t) = A; sin (jrpt/2) / sin (jrp/2), j = 1, 2, and rp = 7r/ (2/3 + 1). Let us consider 
some behaviors of (21) for (3 = f3 r ^ g) where P r>g (g — 1, 2) are the resonant values of 
(3 determined in Sec. Ill: (a) For /3 rjl = r/l — 1/2 mod(l), r = 0, 1,..., lo — 1, one 
has \ki(/3 rj i, t)\ = \k2{(3 r ,i, t)\ = kt. Now, the Bessel function J m (x) is relatively small 
for \m\ > \x\ and, for \x\ < 1, Jo(x) = 0(1) [13]. We then see that on the time scale 
t < T 7 = (Ar/I -1 one can approximate (21) at p = n + /5 r>1 by P(p, t) Jn(kt). Thus, the 
most prominent g = 1 peaks of the MPD for t < T 7 are those with |n| < As 7 — > 
(T 7 — > 00), the expression (21) reduces exactly to P(p, t) = J 2 t)] (see note [14]). (b) 

For (3 r . 2 = r/(2Z )-l/2 mod(l), r = 1, 3, . . . , 2Z -1, one has |fci(/3 r , 2 , t)| = /c or for i odd 
or even, respectively, and \k 2 (P r ,2, t)\ = kt. Then, if \k\ < 1, J n -2m [ki((3 r ,2, t)] is relatively 
small for \n — 2m\ > 1 and one can approximate (21) at p — n + (3 r;2 by P(p, t) m J? n / 2 ](k^t), 



CO 



2 




(21) 



m=— 00 
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where [n/2] is the integer part of n/2. The new (g = 2) peaks thus start to emerge in a 
significant way when t > |/c7| _1 and, for \n\ < 2 j /cyi j , their magnitude should be comparable 
to that of the g — 1 ones. All these behaviors are illustrated in Fig. 1 for r = 2n (l = 1), 
t = 100, k = 0.1, 7 = [Fig. 1(a)] and 7 = 0.2 [Fig. 1(b), with 2\kyt\ = 4]. We have 
checked numerically for many values of (t, k, 7) that the MPD is robust under sufficiently 
small perturbations of r, r = 2n + e, at least within the domains of p where the principal 
peaks above are found. As an example, compare Fig. 1(c) with Fig. 1(b). 

VI. SUMMARY AND CONCLUSIONS 

In conclusion, the results in this paper should provide new insights into the nature of the 
spectra and quantum dynamics of the KP under general QR conditions and for arbitrary 
potentials. As a direct consequence of translational invariance, any rational value of the 
Bloch quasimomentum (5 may be resonant (exhibits QR). At fixed r/{2n) = lo/qo (Iq and go 
are coprime integers), a rational f3 is characterized by an integer pair (r, g) which determines 
(3 through Eq. (7), with I = gl and q = gq assuming their minimal values. The QE 
spectrum consists of q — gq bands which should be, typically, not all flat, implying QR. By 
slightly varying (5 on the rationals, the corresponding value of g changes erratically, leading 
to similar changes in the QE spectrum. The important case of the main QRs (go = 1) is 
precisely described by the linear KR and is thus exactly solvable. In this case, all rational 
values of (5 are indeed resonant for generic potentials containing all the harmonics. In general, 
the resonant values of (3 correspond to peaks in the momentum probability distribution (see 
Fig. 1) and appear explicitly in the formulas (11) and (20) for the diffusion coefficients. 
Formula (20) implies a new phenomenon for multi-harmonic potentials [thus excluding the 
standard case of V(x) = cos(x)]\ A sensitive dependence of the diffusion coefficient on the 
specific localization features of the initial wave packet in momentum space. Assuming the 
robustness of our results under small variations of r, which was numerically verified for 
a two-harmonic potential, it may be possible to observe this phenomenon in experimental 
realizations of the system. While high-order (g > 1) QRs appear to be presently beyond 
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experimental observation, an interesting question is whether cases of such QRs are exactly 
solvable, at least to some extent, for the (5- dependent QE spectra and quantum dynamics. 
We do not have yet a definite answer to this question. We hope that we shall be able to 
make some progress in this direction in future works. 
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APPENDIX 

We derive here the asymptotic behavior (19). The dominant contributions to the in- 
tegral (18) come from small /3-intervals around the zeros of the denominator of the inte- 
grand. We show below that the behavior (19) is completely due to the simultaneous zeros of 
sin(mra/2) and sin(m / r /3 /2) (m, m' ^ 0), Tg = irl (2/3 + 1). The zeros of, say, sin(mr /3 /2) are 
P = Pr,m = r /(\ m \lo) — 1/2 mod(l), r = 0, . . . , \m\l — 1. Writing m = gm and w! = gm' , 
where mo and m' are coprime integers and g = g(m, m') is the greatest common factor 
of (\m\, \m'\ ), it is easy to see that there are precisely gl simultaneous zeros, given by 
P = Pr, g = r/{gl ) - 1/2 mod(l), r = 0, . . . , gl - 1. 

Thus, let j3 = f3 r>m be a zero of sin(mr /3 /2) which is not a zero of sin(m'rg/2), \m'\ ^ \m\, 
and consider a /^-interval I e = \j3 — e, /3 + e] sufficiently small that no zero of sin(m'rg/2) 
lies within it. We show that the contribution of I e to (18) is finite in the limit of t — > oo. 
Let \m\irl e 1, so that sin(mrg/2) (— l) r mTil (f3 — (3) in I e . We assume that the 
correlation function (16) can be expanded as a Taylor series around (5 — (5: Cp{m) = 
Cp{m) +C'p(m)(f3 — ft) H — •. The dominant contribution of I t to (18) is then approximately 
given by 
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B(I n „ (-l) r Cg(m-mQ sin(mr^/2) sin(m V/2) i(ro >_ ro)(t+1)T<)/2 , . 

1 ^ ~ rnvr/o ^ /3 - sin(m'r /3 /2) 1 ; 



We introduce the variable z = irl ((3 — f3)t and define z = irl ((3 + l/2)t. For f>e \ one 

can see that B(I € ) in Eq. (22) is well approximated by 

'-lyft+VCjjjm- m 1 ) f°° dz : _ 
mnlo sin(m'Tg/ 

The integral in Eq. (23) can be calculated exactly using simple trigonometry and formulas 
(3.741.2) and (3.763.2) in Ref. [13]; its value is finite for all z (or t). 



w ( lY^C^m m') /- dz g . n + eiK _ m)M _ 3 

mnlo sm(m / T«/2) J.^ z 



Let us therefore consider a simultaneous zero /3 = /? r)S of sin(mr J g/2) and sin(mVg/2), 
denoting it again by /3. The interval I e is defined as above with 7il e <C min (|m| _1 , |m'| _1 ). 
The contribution of I t to (18) is, approximately, 

* f - m>i (w- m , (t+1 ,^, (24) 

mm 7r t ^ _ 

where m and m(, are the integers defined above. Expanding again Cp(m — m') as a Taylor 
series around (5 — /3, we see that the first-order term C'p(m—m')(f3—j3) in this expansion gives 
a contribution to (24) which is similar to (22) and is thus finite for all t. The contributions 
of higher-order terms are, obviously, also finite for all t. After simple algebra, one can easily 
verify that the contribution B^(I e ) of the zero-order term Cp(m — m!) is well approximated, 
for t > e _1 , by 

(J-(ui m')t f°° 

B {0) (/ e ) « r— / dz z- 2 sm(mz) sin(m'z) cos \(m - m>] , (25) 

mm'ixl J x 

where the variable z was defined above. By expressing the product sin(m',2) cos [(m — m')z] 
as a sum, the integral in Eq. (25) can be calculated exactly using formula (3.741.3) in Ref. 
[13]. The final result, with (3 = j3 r ^ g , is 

B {0) (I e )^ t — (\m\ + \m'\ - \m-m'\)Ca (m-rri). (26) 

2mm'/ 

Summing (26) over all the gl simultaneous zeros f3 = f3 r>g of sin(mr /3 /2) and sin(m'rg/2), 
we obtain the asymptotic behavior (19). 
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FIGURES 



FIG. I. Momentum probability distribution (21) for t = 100, k = 0.1, and (a) r = 2n, 
7 = 0; (b) r = 2vr, 7 = 0.2; (c) r = 2tt + e [e = 7r(V5 - 1)/1200 w 0.0032], 7 = 0.2. All the 
resonant peaks for 7 = [in (a)] are at p = n + 1/2. For 7 = 0.2 [in (b) and (c)], new peaks 
appear at integer p = n. 
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